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Abstract 

We discuss the capabilities of a future neutrino factory based on intense 
neutrino beams from a muon storage ring to explore the non-standard neu- 
trino matter interactions, which are assumed to be sub-leading effects in the 
standard mass induced neutrino oscillations. The conjunction of these two 
mechanisms will magnify fake CP violating effect in the presence of matter 
which is not coming from the CP phase in the neutrino mixing matrix. We 
show that such fake CP violation can be observed in neutrino factory exper- 
iments by measuring the difference between the neutrino and anti-neutrino 
probabilities. In order to perform such test, we consider three neutrino flavors, 
admitting the mixing parameters in the range consistent with the oscillation 
solution to the atmospheric and the solar neutrino problems, as well as with 
the constraints imposed by the reactor neutrino data. We show that with 
a 10 kt detector with 5 years of operation, a stored muon energy En > 20 
GeV, 2 x 10 20 muon decays per year, and a baseline L ~ 732 km, such a 
neutrino facility can probe the non-standard flavor changing neutrino inter- 
actions down to the level of (10 -3 — 10~ 2 ) Gf, in both — > v T jv Vti — ► v T and 
z/ e — ► v T /v e — ► v T modes. 
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I. INTRODUCTION 



The atmospheric neutrino experiments [U present today compelling evidence in favor 
of disappearance. The solar neutrino experiments as well, are strongly indicating v e 
disappearance. The Los Alamos Liquid Scintillator Neutrino Detector (LSND) || has also 
reported results consistent with neutrino flavor conversion. 

The most plausible mechanism of such neutrino flavor disappearance or conversion is 
the oscillation of neutrino due to quantum interference driven by neutrino mass squared 
differences Am 2 ||||, which can happen when neutrino flavor eigenstates are supposed to 
be coherent superpositions of neutrino mass eigenstates |J. Since the impressive results from 
the LEP experiments [§|,|7J|, supported recently by the direct observation of v T by DONUT [g, 
established that there are at least three active neutrino flavors, it is unavoidable to try to 
understand such mass induced neutrino oscillations in a full three generation scenario. 

Under the assumption of such mass induced neutrino oscillation in three flavor frame- 
work, atmospheric neutrino data || can be explained by pure — > v T oscillations in vac- 
uum |10[ , though the possibility of having contributions from non- negligible z/ M — > v e oscil- 
lations is still not discarded 0] , even after taking into account the constraints coming from 
the CHOOZ reactor experiment [12|]. On the other hand, solutions to the solar neutrino 
problem can be provided either by the matter enhanced resonant neutrino conversion, the 
MSW effect fli3||U| , or by vacuum oscillations with a typical wavelength of the order of the 
Sun-Earth distance |15| with the pure two generation mixing or with small additional par- 
ticipation of a third neutrino on top of the two generation mixing. In this work, we do not 
take into account signals of the LSND experiment, since simultaneous explanations of the 
LSND data together with atmospheric as well as solar neutrinos data require the presence 
of a fourth sterile (electroweak singlet) neutrino, which is beyond the scope of this paper. 

Several non-standard (or exotic) explanations of atmospheric as well as solar neutrino 
data, which do not necessarily invoke neutrino mass and/or mixing, also have been sug- 
gested. For atmospheric neutrinos, conversion or disappearance mechanisms induced 
by neutrino decay [TB|, flavor changing neutrino interactions |ff7| and quantum decoher- 
ence |18| have been proposed to explain these data. Among them, it was found that the 



flavor changing solution can not explain well the SK upward going muon data |n| and also 
it is disfavored by the results coming from K2K f2~0|| . For solar neutrinos, there are also 
several proposals to explain the data by invoking non-standard neutrino properties such as 
magnetic moment ETJ, non-standard neutrino interactions [55 or even by a tiny violation 
of the equivalence principle [^3|] . 

While some of these non-standard solutions can explain quite well the data, it is gener- 
ally believed that the most plausible solutions are provided by the standard mass induced 
oscillation mechanism, because mass and mixing are the simplest extension of the standard 
model and moreover, it is the unique mechanism which can explain well both atmospheric 
and solar neutrino observations at the same time, without additional neutrino properties. 

However, even in this case, we consider that it is quite possible that some non-standard 
neutrino properties could induce some sub-leading effect in addition to the the standard mass 
induced neutrino oscillation without causing any inconsistency with the present observations. 
We are particularly interested in such sub-leading effects induced by flavor changing neutrino 
interactions. These type of interactions were first suggested in Ref. 
as a possible mechanism to solve the solar neutrino problem in Ref. 



2~3}| , and then discussed 



25.26 



Indeed, many models of neutrino mass p7H29| are also a natural source of non-standard 



2 



neutrino matter interactions (NSNI) which can be, phenomenologically, classified into two 
types : flavor changing neutrino interactions (FCNI) and flavor diagonal neutrino interac- 
tions (FDNI), as we will discuss in detail in our work. 

Therefore, we believe that even if mass induced oscillations really take place in nature 
and such NSNI are not playing any relevant roles in explaining the observed neutrino data, it 
could be important to study NSNI interactions as secondary contributions, in order to gain 
some handle on new physics beyond the electroweak standard model. The crucial point of 
our idea is to use the fact that the simultaneous presence of neutrino masses and NSNI can 
enhance the difference in the conversion probabilities between neutrino and anti-neutrino 
channels in matter causing fake CP violation |30| . 

Our aim here is to test such interactions in a future neutrino factory. A neutrino fac- 
tory |3~lH , which could be a milestone for a future muon-collider project, presents a possibility 
to go beyond the goal of measuring with great precision the oscillation parameters. It will 
offer us the chance to measure separately the oscillating probabilities in CP-conjugate chan- 



nels, such as v.. 



muon storage ring [31 



and v 



with high intensity neutrino beams produced by a 



after neutrinos travel through several hundreds of kilometers in the 
Earth mantle before reaching the detector. 

There are already quite a few very interesting studies that have shown the prospects of 
neutrino factories to measure the mixing parameters and obtaining the pattern of neutrino 
to test matter effects 



IM3- 



32,33], CP violation in the leptonic sector 



masses 

CPT-odd interactions |36| and distinguish the mass induced solution for the atmospheric 
neutrino anomaly from the decoherence one |J7[] , or possible neutrino oscillation induced by 
violation of the equivalence principle [38]. Signals of R— parity violating supersymmetric 
interactions which can be mistaken for neutrino oscillations in neutrino factory experiments 
were also investigated in Ref. 



39 



where it is found that these novel contributions basically 
do not affect the r event rate for baselines greater than 200 km, except if the new couplings 
are close to their perturbative limit. 

Here, we investigate the capacity of neutrino factory experiments to probe neutrino 
oscillations beyond the standard mechanism assuming mixing parameters compatible with 
atmospheric, solar and reactor neutrino experiments and by so doing to discover new physics 
or put restrictive limits on NSNI models. Our results will be presented as a function of the 
distance between source and detector for two particular muon energies. 

In Sec. |TJ, we define the formalism that will be used throughout this paper. In Sec. [JTI|, 
we describe the various experimental set-ups we will consider, as well as the details of our 



experimental simulations. In Sec. [TV| , we show our predictions for the 



and v e 



channels with and without CP violation, 
present our conclusions. 



Finally in Sec. [V|, we discuss our results and 



II. MODIFIED OSCILLATION FORMALISM WITH NSNI 

There are currently several experimental bounds on NSNI coming from the non- 
observation of lepton flavor violating process |J] , from violation of lepton flavor universality, 
or from neutrino scattering data. Model independent analysis of these experimental con- 
straints derived from SU{2)l related interactions, that can be found in Refs. [PC| , ^1 , in 
general, gives stronger bound for v e — channel, typically set the maximal strength of 
NSNI normalized by the electro weak interactions strength, i.e. Gf, to be at the level much 
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smaller than one-percent whereas the bounds for — v T and u e — v T channels are substan- 
tially weaker typically at the level of ~ 10~ 2 — 10 _1 in unit of Gp. Therefore, in this work, 



we will only consider the two channels za< 



and their anti-neutrino ones, which 



are less constrained from the laboratory experiments. 

In the usual three generation neutrino oscillation framework, we can define the corre- 
spondence between neutrino mass eigenstates and neutrino interaction eigenstates by, 



(2.1) 



where Ui (i = 1, 2, 3) are the neutrino mass eigenstates, and v a (a = e, /i, r) are the neutrino 
flavor eigenstates with the standard parameterization for Maki-Nakagawa-Sakata || mixing 
matrix U given by, 
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(2.2) 



•S12S23 — c 12 c 23 s 13 e — C12S23 — Si2C23Si3e C23C13 

where c^- and Sij denote the cosine and the sine of dij and 5 is the CP violating phase. 



A. NSNI in v„ 



channel 



If we consider the possibility of NSNI only in the 
tonian in matter has the form, 



v T channel, the evolution Hamil- 
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(2.3) 



{mf — m/j)/2E, i,j = 1,2,3, E is the neutrino energy, V/(r) = 
is the flavor-changing + f — *■ v T + f forward scattering 



amplitude with the interacting fermion / (electron, d or u quark) and Vf(r)e'f is the differ- 
ence between the flavor diagonal v^ — f and u T — f elastic forward scattering amplitudes, 
with rif (r) being the number density of the fermions which induce such processes. In the 
evolution Eq. ( |2.3| ), e^ T and e'f are the phenomenological parameters which characterize 
the strength of FCNI and FDNI, respectively. The fermion number density n./(r) can be 
written in terms of the matter density p as n/(r) = p(r)Yf, where Yf is the fraction of the 
fermion / per nucleon, ~ 1/2 for electrons and ~ 3/2 for u or d quarks. In all cases, electron 
neutrinos will coherently scatter off the electrons present in matter through the standard 
electroweak charged currents which introduce non trivial contributions to the neutrino evo- 
lution equations. These contributions are taken into account by the diagonal term V e (r) in 
Eq. (p|). 

A similar form applies to the evolution of anti-neutrinos except that V e j(r) —> —V G j(r) 
and U — > U*. We assume the density profile of the Earth to be the one given by the 
Preliminary Reference Earth Model and solve these evolution equations numerically to 
compute the oscillation probabilities Pivp, — > 



and P(v,, 
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B. NSNI in v e — * v T channel 



Now, if we consider NSNI only in the v e 
matter will have the form, 



d 



v T channel, the evolution Hamiltonian in 
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(2.4) 



Again for anti-neutrinos we have to make the following substitution V e j{r) — > —V e j(r) and 
U — > U* . We obtain the oscillation probabilities P[y e — > u T ) and P{p e — * u T ) by numerically 
solving Eq. ( |2.4| ) for neutrinos and the corresponding modified one for anti-neutrinos. 

In the following, we will drop the / superscript of e, e' and consider interactions only 
with either u or d quarks, since limits on interactions with electrons can be obtained by a 
simple re-scale of our plots. If NSNI are caused by electrons, e and e' parameters must be 
increased by a factor 3 to get the same effect presented in this paper. 

C. Useful Two Generation Formulae 



It is instructive to write down the analytical expression for the conversion probability in 
two generation considering constant matter density. This will help to understand the full- 
fledged three generation numerical results. The probability of v a — > u T , a = e,fi, conversion 
can be written as 



P(v n 



'A n 2 

I ( -sm26 + e aT V f 



4 ( S i n 26 + e aT V f ) 2 + (A cos 26 + e' ar V f - 6 ea V e 



L 

x sm" < — 
2 



sin 26 + e ar V f f + (A cos 26 + e' ar V f - 5 ea V e 



(2.5) 



where A and sin26 l will be chosen depending on the oscillation mode. The probability of 
u a — > v T conversion is analogous to the above except that V e j — > —V e j. 

One can easily see that the probability is much more sensitive to e aT than to e' aT , since 
e aT appears in the numerator and in the denominator of the amplitude in Eq. (|2.5| ), while 
e' QT appears only in the denominator. Moreover, the terms involving e aT are increased by a 
factor 4 with respect to the terms involving e' ar . These features become more prominent if 
the distance is not so large and the contribution from FCNI are small as we will see below. 

Let us take some typical values for the mixing parameters and the neutrino energy we 
are going to consider in this work, Am 2 ~ few x 10~ 3 eV 2 , E = few xlO GeV and not 
very large distance, L < 1000 km. In this case, taking the NSNI effects to be small, i.e. , 
e,e'<l, the probability in Eq. ( |2.5|) can be approximated, up to the first order in the NSNI 
parameters, as follows, 



P(Vn 



1 



A sin 26 (A sin 2^ + 4 e aT V f ) L 2 



(2.6) 



Note that there is no dependence on the e' ar parameter in the above probability. For the 
corresponding anti-neutrino channel, Vf must be replaced by —Vf in Eq. ( |2.6| ), which lead 
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to a different probability. It is important to point out that, as we can see from Eq. ( j2.5|) , 
only the presence of NSNI term alone without neutrino mass (A — > 0) can not cause the 
neutrino and anti-neutrino probabilities to differ. 

In order to extract the information on NSNI, without knowing the very precise values 
of mixing parameters, we must compare in some way, the probability for neutrinos and 
anti-neutrinos. In this work, we consider the ratio of the expected number of events for 
anti-neutrinos and neutrinos, which can be approximately inferred, for small values of e aT 
taking into account the difference in the detection cross sections [o Dt and cr„ r ), from the 
simple expression, 

R(p) = ^ T P(i? Q ~+ Vr) „ 1 _ 4 fgrVf (2 7) 

o VT P{v a -»• u T ) 2 A sin 29 + 4 e aT V f ' K ' J 



From the expressions in Eqs. (|2.6| ) and ( |2.7| ), it is clear that the oscillation probabilities are 
much more sensitive to the FCNI parameter, e aT , than to the FDNI parameter, e' aT , and 
moreover, in the ratio, the FCNI effect does not depend on the distance as long as it is not 
very large, i.e. L < 1000 km, this will be confirmed in the following sections. 

In this work, we assume that Am^ is positive and do not consider the negative case 
(inverted hierarchy) because such case will covered as we consider both the positive and 
negative sign of e ar . /,From Eq. (|2.7|) we can see that the same effect is obtained when 



(Amjj, e aT ) is replaced by (-Am^, — e aT ), which is valid as long as L < 1000 km. 



III. EXPERIMENTAL SET-UP AND OBSERVABLES 

Many authors [31] have emphasized the advantages of using the long straight section of 



a high intensity muon storage ring to make a neutrino factory. The muons (anti-muons) 
accelerated to an energy constitute a pure source of both (z/ M ) and v e iy e ) through their 
decay \l~ — > e'VeV^ — > e + v e v^) with well known initial flux and energy distribution. 
In this context it is quite suitable to perform very precise measurements of the probability 
of oscillation in CP-conjugated channels, such as — > v T and — > v T or v e — > v T and 

There are many propositions for these type of Neutrino Factories, with values for the 
energy of the stored muon (anti-muon), E^, going from 10 GeV to 250 GeV, for baselines, 
L, ranging from 730 km to 10 000 km, in any case the neutrino beam penetrating a fair 
bit of the Earth's crust. Since we do not know which will be the final configuration we will 
do our estimations for several possible configurations. We do our calculations taking into 
account all the available relevant experimental information to compute real experimental 
observables. 

Here we will explore the neutrino factory as an appearance v T (u T ) experiment. We define 
the following observables of interest: 

iV i/ T v e / iT ' Lit) 



and 



^-4) = §tHi' (3 ^ 2) 
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where (R er ) is the ratio between the total number of detectable z/ T -events, N£ (N^), 
when the neutrino beam is made of (fi~) decays, in the — > v T {p e — > v T ) channel, over 
the total number of detectable ^ T -events, N~ t (iV+), when the beam is made of 
decays, in the — > v T [y e — > v T ) channel. These numbers can be calculated as 

10 9 Na E 3 r E » 

N-ie^e'^) = ^_M— ^-f \ h{E)P{v, -> u T )(E)dE, (3.3) 

1 D 9 AAa E 3 r E u, - 

*£(^,<J = n^-M— r^T? / HE)P(u, - u T )(E)dE, (3.4) 

m^TT L, JE th 



11) iv A Hi /"^M 

N+ T (e eT ,e' eT ) = n M -M^^-f / ^(E)P(i/ e - u T )(E)dE, (3.5) 

Tfl^TT Li J E th 

10 9 Na E 3 r E » 

N^(e eT ,e' eT ) = n^Af— ^-f / g(E)P(u e - z/ T )(£)ciE, (3.6) 

m u' K ^ JE ^ 



with 



h(E) = 


E\ 


-2# 




h(E) = 


4r 


-2—; 




9(E) = 


E 2 
12— (1 

Et 


E 
~% 




9(E) = 


-fa 


E, 
E/ 





(3.7) 

(3.8) 

(3.9) 

(3.10) 
(3.11) 

where is the muon source energy, M is the detector mass in ktons, n M - and the 
number of useful /i~ and /i + decays, respectively, 10 9 ./Va is the number of nucleons in a kton 
and m M is the mass of the muon. The functions h(E) and h(E) (g(E) and g(E)), contain 
the and [y e and z/ e ) energy spectrum normalized to 1, the charged current interaction 
cross section for v T and v T |42]| , and the u T , v T detection efficiencies r] T ,ri f , respectively. 



We have used E t h = 4 GeV and assumed, for simplicity, that once this cut is applied the 
tau neutrino and tau anti-neutrino, can be observed with the same efficiency r] T = r} f = 0.33 
which seems to be achievable by means of exploiting one-prong and three-prong r decay 
topologies along with displaced vertex or kinks resulting from r-lepton decays f43|. We have 
neglected the finite detector resolution following Ref. |3J2] . 



In our study we will consider the following general characteristics for the neutrino factory : 
E^ = 20 GeV and 50 GeV, 2 x 10 20 muon decays, a 10 kton detector and 5 years of data 
taking. We would like to point out that the observable R^ T (R er ) can amplify the differences 
between iV~(iV+) and N^ t (N^_), independent of the absolute number of events. 
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IV. PERSPECTIVE OF FUTURE NEUTRINO FACTORIES 



In order to estimate the maximal limit on e aT and e' aT , with a = e, fi that can be achieved 
by a neutrino factory we define the following functions 



,2/, J \ _ nfAT± 




2[<(6c 



N±(e ar ,e' aT ) 



(4.1) 



and 




X 2 (e aT ,e' ar ) = 2{Nl(e ar ,e' a 

+ 2Nl(e aT = 0,e' QT = 0) In I ^ UtK "JZ ( ^ ] , ( 1.2) 

so that for fixed values of the oscillation parameters, one can compute the number of standard 
deviations of separation between pure mass induced oscillation and mass induced oscillation 
plus NSNI contribution as a function of e aT and e' aT , as n a = J (x 2 + X 2 )/2. 

We will consider in our analyses the following range for the oscillation parameters : 
1 x 10~ 3 eV 2 < lAmlJ < 7 x 10" 3 eV 2 , \Am 2 21 \ < 1 x 10" 3 eV 2 , 0.8 < sin 2 26 23 < 1.0, 
0.7 < sin 2 20i2 < 1.0 and sin 2 2#i 3 < 0.3; with < 9 12 , 9 23 , 9 13 < k/2 and < 5 < vr/2. To 
the extent that the contributions from the sub- leading solar [Am 2 .^ << jAm^l scale are 
small, our results apply approximately to other solar scenarios. 



A. NSNI in — > v T channel 



The interesting feature of this channel is not only that for the standard oscillation mech- 



anism, with 5 = 0, the fake CP violation induced by matter effect is very small fl44] , under 
our assumptions about the mixing parameters, but also that even if we consider the CP 
violation phase to be maximal (S = vr/2), this quantity remains quite negligible. This is 
due to the fact that sin 2 2#i3 is highly constrained by CHOOZ and the atmospheric neutrino 
data. On the other hand, as we have seen in Sec. |I| the inclusion of an extra contributions 
from NSNI in the neutrino evolution Hamiltonian will enhance the fake CP violation in 
matter. 

To illustrate the size of the effect due to NSNI in this channel we plot in Fig. [T] and [2] the 
ratio R^ T as a function of L for = 20 and 50 GeV respectively, for e^ r = and e MT = 0.001 
and 0.01. These plots were done for the best fit values of the oscillation parameters for the 



combined analysis of atmospheric, solar and reactor data according to Ref. |4"5 |, assuming 
the large mixing angle solution to the solar neutrino anomaly, i.e. , Am^ = 3.1 x 10~ 3 eV 2 , 
Am 2 x = 3.3 x 10" 5 eV 2 , sin 2 28 23 = 0.97, sin 2 2# 12 = 0.78 and sin 2 28 n = 0.02, with 5 = 0. 

For no NSNI, the ratio is nearly constant around 0.5, it is essentially dominated by the 
ratio between v T and v T cross sections since the standard matter effect is negligible {Piy^ — > 
u T ) w P{Vfx — > u T )) . As FCNI are introduced by increasing the value of e^ T , we observe 
the effect of the new interactions pop in, becoming quite appreciable for e Mr ~ 10~ 2 and 
stronger as we increase in energy. This dependence in energy can be easily understood using 
the two generation formula given in Eqs. ( p.5|) -( ^7|) , setting A = A 32 and sin 29 = sin2# 2 3- 
Although we are here working in a full three generation scheme, this approach will be good 
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enough to interpret the general behavior of the probabilities. This is because these two 
parameters are the leading parameters for conversion in this channel in the three generation 
framework. ^From Eqs. (|2.6|) and ( 2/7) , we see that the effect of fake CP violation become 
larger as energy grows and can reach smaller values of e Mr . 

Note that the impact of using the Earth's density profile can be appreciated in the small 
distortions of the curves for L > 1000 km, which is specially sizable in Fig. |2], for e MT =0.01, 
which also can be explained in first approximation by the two generation formula. 

For our statistical analyses of the limits that can be achieved by a neutrino factory, we 
have chosen to work with L = 732 km. Two reasons support this choice : the FCNI effect is 
stronger and independent of L up to 1000 km (which can be understood from Eq. (|2.7|) in 
Sec. [II C|) and this baseline is compatible with the CERN-Gran Sasso or Fermilab-Soudan 
distance. Besides this at 732 km the neutrino flux is still quite big. 

In Figs. and f|, we show the region of sensitivity in \e'\ x je^l plane for e' , > 
computed using the prescription given in Sec. |V| and requiring a n a = 3, 5 separation. We 
plot the maximal limit for two extreme cases, one for Am^ = 1 x 10 -3 eV 2 and the other 
for Am 2 2 = 7 x 10 -3 eV 2 . To give an idea of the corresponding number of events expected 
at the limiting points, we give in Tab. [I] these numbers for e^ r = 0. We have checked that 
these results are not modified by changes in the values of the other oscillation parameters or 
by the adopted sign of Am§ 2 or of Am 21 . In particular, P{v^ — > u T ) and the P{v^ — > u r ) are 
essentially independent of 5 and the value of Am 21 . We also have done the same calculation 
in the case of either e^ T < or e AtT < or both, but we do not show these curves here since 
they are almost identical to the ones in Figs. ^ and |j. 

For En = 20 GeV we see that it is possible to probe, at 3 o level, |e MT | > 6 x 10~ 3 
and le^l > 1.7 (0.6), for E^ = 50 GeV, |e MT | > 1.6 x 10~ 3 and |e^ T | > 1.4 (0.4), for 
Am 2 2 = 1 x 10 -3 eV 2 (7 x 10~ 3 eV 2 ). As explained in Sec. [TV] we did not expect to be very 
sensitive to e^ r , also the limits are better at 50 GeV due to the flux increase with E^. 

Finally we note that, although the oscillation probabilities in this channel can get quite 
large, the difference in the neutrino and anti-neutrino probabilities remains small, unless 
\ e fir\ 1- This means that the constraints relies mainly on a large event statistics. 



B. NSNI in u e — > v T channel 

We plot in Figs. [5] and |6] the ratio R eT as a function of L for E^ = 20 GeV and 50 GeV, 
respectively, for e' eT = and e er = 0.0001,0.001 and 0.01. We can compare these figures 
with Figs. [I] and |2] for and make two observation. First, for the case without NSNI 
we have a strong deviation in the channel v e — > v T from the value 0.5 after L > 1000 km, 
this is because matter effects (V e ) are specially important in this mode, in contrast with 
Vfj, — > v T where they are small. Second, when we add the non-standard contributions the 
effect in v e v T channel is again more relevant than for — > v T . This can be understood 
since at first approximation we can think that the modified three generation probabilities 



P{y e — > v T ) and P{p e — > z/ T ), as it is the case for the standard ones p3l , must be proportional 
to sin 2 29™ 3 in matter, only in our case we also have to consider the contributions of the non- 
standard interactions. This permits a qualitative understanding of our results in terms of 
the two-generation formula (Eq. ( |2.5[ )) replacing A by A 31 and sin 28 by sin 2#i3 , which are 
the most relevant parameters in this mode. ^From Eqs. ( |2.6| ) and ( |2.7| ), we see that the 
relative effect of fake CP violation become larger compared to the z/„ — > v T channel because 
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of the smaller value of the relevant mixing angle 6i 3 for this channel. 

In Figs. |?|-|I0|, we show the minimal values of e eT that can be probed at = 20 and 50 
GeV. The first two figures correspond to the cases e er , e' er > and e er , e' eT < at 20 GeV, 
the last two to the cases e er , e' eT > and e er , e' eT < at 50 GeV. We also have looked at the 
possibility of having e er > 0, e' eT < and e er < 0, e' eT > 0, the limits in the first case are 
almost identical to the ones shown in Figs. |7| and ^ while in the latter case very similar to 
Figs. | and 0. 

These plots were calculated again assuming a comfortable separation of n a = 3, 5 in 
the x 2 analysis described in Sec. [IV] and L = 732 km as for the — ► v T mode. For 
each plot we have two different assumptions (a) 5 = and (b) 5 = n/2, in both we have 
considered the uncertainty on Am^. We see clearly the CP violation phase 5 is significant 
in this channel. We do not show the dependence of these limits on the values of the other 
oscillation parameters or on the sign of the mass squared differences, since these would 
practically unaffect the plots. In the Tab. |, one can find the corresponding number of 
events when e eT > and e' eT = for n a = 3. 

For = 20 GeV, we see that it is possible to probe, at the 3 a level, for any value of 
e'er; £er ^ 7 x 10~ 3 and e er < — 7 x 10~ 3 for 5 = 0. For 5 = tt/2 one can test e er > 1.2 x 10~ 2 
( 1.5 x 10~ 2 ) and e er < —1.4 x 10~ 2 ( —5.1 x 10 -2 ), these limits correspond to Am 2 2 = 
1 x 10~ 3 eV 2 (7 x 10~ 3 eV 2 ). 

For = 50 GeV, it is possible to test e er > 1.6 x 10~ 3 and e er < —1.6 x 10~ 3 for 5 = 
and e eT > 4.0 x 10~ 3 ( 6.3 x 10~ 3 ) and e er < -3.6 x 10~ 3 ( -1.4 x 10" 2 ) for 5 = tt/2, again 
these limits correspond to Am 2 2 = 1 x 10~ 3 eV 2 (7 x 10~ 3 eV 2 ). In this channel it will be 
not possible to set a stiff limit on e' eT due to low statistics. 

In spite of the low statistic of events in this mode, the difference in the probabilities of 
oscillation for neutrinos and anti-neutrinos in matter, as one can see from Figs. || and ||, 
makes it possible to achieve a sensitivity to FCNI as low as the previous mode. 



V. FINAL CONCLUSIONS 

We have studied the effect of non-standard neutrino matter interactions which are as- 
sumed to be sub-dominant in the standard mass induced neutrino oscillation in a future 
neutrino factory experiments. When non-standard matter interactions is added to the mass 
oscillation mechanism, the difference between neutrino and anti-neutrino probabilities can 
be enhanced, causing fake CP violating effect not coming from the intrinsic CP phase. 

Based on this interesting feature, we have obtained, in v e —>■ v T as well as — > v T 
oscillation modes, possible limits on NSNI which can be accessible by the future neutrino 
factories, under the frame work of three neutrino mixing scheme, assuming the mixing 
parameters compatibles with the solar and atmospheric solution as well as reactor neutrino 
oscillation experiments. 

We have shown that these facilities, with 5 years of operation and a 10 kt detector, are 
able to test FCNI down to a level of (10~ 3 — 10~ 2 ) Gp in either z/ M — * v T and v e v T modes 
whereas sensitivity to FDNI is much worse only to a level of 0(1) Gp for — > v T . These 
limits on FCNI will certainly put much more stringent limits on the strength of lepton flavor 
violating new interactions, parameterized here by e, than the present bound can be found 
in the literature. • 



In particular, we conclude that if En = 20 GeV and L = 732 km, neutrino factory 
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oscillation experiments will be able to set limits on the FCNI parameters such as |e^_ d | < 
0.006, |e£ T | < 0.018, -0.006 < e u J < 0.007, -0.018 < e% < 0.021, and if = 50 GeV these 
limits can be < 0.002, |e^ r | < 0.006, -0.014 < e u J < 0.007, -0.042 < t e eT < 0.021. 
These bounds could be regarded as robust and model independent and can be used to 
constraint new physics in the electroweak sector. Let us also remark that, for both — > v T 
and v e — > v T channels, our bounds are mainly affected by the lAmj^l and do not practically 
depend on [Am^l which is relevant for the solutions to the solar neutrino problem, and 
there is some effect from the CP violation phase, 5, for the v e — > v T channel. 

The constraints on e^ r that could be obtained are e'^ d < 2 and e^ r < 0.7. As we had 
expected these limits are much weaker than the corresponding ones in e^, due to the fact 
that the modified probability (mass+NSNI) is slightly dependent in e^ T (see discussion on 
sec II. C) for the distance consider. Also these constraints are clearly much less stringent 
than the ones in Refs. [4C,^]. Due to low statistics no limit on e' eT can be established. 



We have not present the analysis for e eAt since it will be very hard to lower the bounds 
founded in the literature 0(1O~ 5 ) by a neutrino factory. 

As we were finishing this work we came across Ref. |46| where the effect of extra real 



as well as fake CP violation due to new physics was studied in the context of oscillation 



experiments at neutrino factories. We note that the work in Ref. \W\ is somehow com- 
plementary to our since there, the effect of new interactions was considered in the process 
of production and detection of neutrinos, while here we only consider their effect in the 
neutrino propagation. 
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TABLES 



TABLE I. Number of v T and v T events expected for two different values L = 732 km, 2 x 10 
\x decays, assuming a 10 kton detector after 5 years of data taking. We give this information 
for Am| 2 = 1 x 10~ 3 eV 2 and 7 x 10~ 3 eV 2 , keeping the other oscillation parameters fixed to 
Am 2 .! = 3.3 x 10~ 5 eV 2 , sin 2 2# 23 = 0.97, sin 2 29 12 = 0.78, sin 2 29 13 = 0.02 and 5 = 0. Here 
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FIGURES 
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FIG. 1. Ratio between the number of expected v T over the number of expected function 
of the baseline L, for a few values of e MT and e^ T = 0. The value of the other oscillation parameters 
used in this plot are : Amj 2 = 3.1 x 1(T 3 eV 2 , Am^ = 3.3 x 10~ 5 eV 2 , sin 2 2023 = 0.97, 
sin 2 26»i2 = 0.78, sin 2 26 13 = 0.02 and 6 = 0. 
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FIG. 2. Same as Fig. | but for = 50 GeV. 
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FIG. 3. Region in the e^ T x e^ T plane that can be probed by a neutrino factory considering 
En = 20 GeV, L = 732 km, 2 x 10 20 muon decays, a 10 kton detector after 5 years of data taking. 
The regions that were determined by demanding a n a = 3, 5 separation, are the ones to the right of 
each line. The value of the other oscillation parameters used in this plot are : Am^ = 3.3 x 10 -5 



eV 2 , sin 2 29 23 = 0.97, sin 2 26 12 = 0.78, sin 2 20 13 = 0.02 and 5 = 0. Here e^, e'^ T > 0. 
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FIG. 4. Same as Fig. | but for En = 50 GeV. 
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FIG. 5. Same as Fig. |l| but for the mode 
5 = vr/2, -vr/2. 
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FIG. 6. Same as Fig. | but for = 50 GeV. 
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FIG. 7. Same as Fig. |3| but for the mode v e — » v T jv e —* z/ r and two distinct values of S. Here 
e er ,eg T > 0. 
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FIG. 8. Same as Fig. |3| but for the mode v e — » v T jv e —* z/ r and two distinct values of <5. Here 
e er ,eg T < 0. 
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FIG. 9. Same as Fig. @ but for £ M = 50 GeV. 
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FIG. 10. Same as Fig. | but for £7„ = 50 GeV. 
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